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FOREWORD 


This report was prepared by North American Aviation, Inc. , Space 
Division, under NASA Contract NAS9-4552, for the National Aeronautics and 
Space Administration, Manned Space Flight Center, Houston, Texas, with 
Dr. F.C. Hung, Program Manager and Mr. P. P. Radkowski, Assistant 
Program Manager. This work was administered under the direction of 
Structural Mechanics Division, MSC, Houston, Texas with Dr. F. Stebbins 
as the technical monitor. 

This report is presented in eleven volumes for convenience in handling 
and distribution. All volumes are unclassified. 

The objective of the study was to develop methods and Fortran IV 
computer programs to determine by the techniques described below, the 
hydro-elastic response of representation of the structure of the Apollo Com- 
mand Module immediately following impact on the water. The development 
of theory, methods and computer programs is presented as Task I Hydro- 
dynamic Pressures, Task II Structural Response and Task HI Hydroelastic 
Response Analysis. 

Under Task I - Computing program to extend flexible sphere using the 
Spencer and Shiftman approach has been developed. Analytical formulation 
by Dr. Li using nonlinear hydrodynamic theory on structural portion is 
formulated. In order to cover a wide range of impact conditions, future 
extensions are necessary in the following items: 

a. Using linear hydrodynamic theory to include horizontal velocity 
and rotation. 

b. Nonlinear hydrodynamic theory to develop computing program on 
spherical portion and to develop nonlinear theory on toroidal and 
conic sections. 

Under Task II - Computing program and User's Manual were developed 
for nonsymmetrical loading on unsymmetrical elastic shells. To fully 
develop the theory and methods to cover realistic Apollo configuration the 
following extensions are recommended: 

a. Modes of vibration and modal analysis. 

b. Extension to nonsymmetric short time impulses. 


iii 


c. Linear buckling and elas to -plastic analysis 

These technical extensions will not only be useful for Apollo and 
future Apollo growth configurations , but they will also be of value to other 
aeronautical and spacecraft programs. 

The hydroelastic response of the flexible shell is obtained by the 
numerical solution of the combined hydrodynamic and shell equations. The 
results obtained herein are compared numerically with those derived by 
neglecting the interaction and applying rigid body pressures to the same 
elastic shell. The numerical results show that for an axially symmetric 
impact of the particular shell studied, the interaction between the shell and 
the fluid produces appreciable differences in the overall acceleration of the 
center of gravity of the shell, and in the distribution of the pressures and 
responses. However the maximum responses are within 15% of those pro- 
duced when the interaction between the fluid and the shell is neglected. A 
brief summary of results is shown in the abstracts of individual volumes. 

The volume number and authors are listed on the following page. 

The contractor's designation for this report is SID 67-498. 
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ABSTRACT 


This volume is a user's manual for a computer pro- 
gram which determines the dynamic response of a shell 
of revolution during a vertical axially symmetric impact 
into an incompressible fluid. The program uses the 
theory developed in Volume 2 of this report where no 
interaction between the fluid and the flexible shell is 
accounted for. The hydrodynamic pressures are deter- 
mined on the basis of a rigid-body theory, and are applied 
to the shell as a forcing function. The results are intended 
for comparison with similar calculations derived from 
Volume 3 and the User's Manual of Volume 6 where the 
hydroelastic interaction is accounted for. 
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1.1 INTRODUCTION 

The computer program for the calculation of the dynamic response 
of shells of revolution during vertical impact into water when no interaction 
is present is written entirely in FORTRAN IV and makes use of the overlay 
feature of that language. The program has been checked out in NAASYS, the 
NAA adaptation of the IBM 7090/7094 IBSYS/IBJOB system; and uses the 
NAASYS library routines shown in the load map, pages 4 to 8, inclusive, 
of Section 1.2. : v ' 

The NAASYS input tape is Unit 5, the output tape is Unit 6. In addition 
to these files, the program uses Units 8, 9, 10, and 11 as scratch tapes, 
and Unit 7 as the overlay tape. NAASYS itself is stored on Unit 1. 

The program is made up of an executive program and eight links, all 
of which are called by the executive program. A brief description of each 
link is shown in Table I below. 


Table 1. Description of Dinks 


Link 

No. 

Name 

Purpose 

0 

Executive 

Reads general data, DA, and controls flow of 
execution of other links 

1 

GEOM 

Reads geometric parameters. Prints all geometric 
input and calculated values 

2 

CDAFIT 

Sets up stiffness parameters 

3 

ACCN 

Computes hydrodynamic pressures on the shell 

4 

DEFLTN 

Calculates the deflections due to the pressures 

5 

PATH 

Controls flow after computation of deflections. 
Computes velocities and accelerations 

6 

INTLDS 

Computes internal loads 

7 

PSUMS 

Outputs all computed quantities 

8 

PIX 

A dummy subroutine for a CRT Plotter 
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2.1 PROGRAM FLOW DESCRIPTION 

An overall flow diagram of the executive program 157DR is shown in 
Figure 1. A listing of the complete program is shown in Section 7. 1. 
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Figure 1 



Flow of Executive Program 157 DR (Sheet 1 of 3) 
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TFI 


Ore. .0 \ / UT*0 



Call DEFLTN 
calculate 
deflections 


Call PATH 

determines program 


Computes 

velocity S accelerations 
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2.2 Program Deck Setup 

As explained in Section 1. 1, the deck is set up in overlay regions. 
Each region is denoted by a $ORIGIN control card. A list of the setup is 
shown below. It includes the control cards and deck names. The order of 
these decks must be kept in the given sequence. 
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•• -M * 

' .V 


\ v.’ *■ v* 

Control Cards 

Subroutines 

$IB JOB 



,:|..$jBFTC 

157 DR 

Main program 


MMPY 

MM Y 

'■H $»FTC 

MADD 


' $ORIGIN 

CHAIN 


$IBFTC 

GMTRY 

GEOM 

$IBFTC 

■ CF3P'.-:--/ ; 

CODIMA 

$ORIGIN 

CHAIN 


$IBFTC 

CDAFIT 

CRVFIT 

49BFTC 

CODS 

CODIMA 

$IBFTC 

ENTP 

ENTERP 

$ORIGIN 

CHAIN 


$IBFTC 

ACCN2 

ACCN 

$IBFTC 

PMAXX 

PMAXL 

$ORIGIN 

CHAIN 


$IBFTC 

157DRI 

DEFLTN 

$IBF T C 

MSUB 

MSU 

$IBFTC 

INVRS 

INV 

$ ORIGIN 

CHAIN 


$IBFTC 

WHERE 

PATH 

$ ORIGIN 

CHAIN 


$IBFTC 

157DR2 

INTL.DS 

$ORIGIN 

CHAIN 


$OBFTC 

FSUMS 

SUMS 

$ORIGIN 

CHAIN, SYSUT2, REW 


$IBFTC 

LNK6 

PIX 

$DATA 
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3.1 RESTART 

In many calculations, it may be desired to restart the program at some 
time t fi without recalculating all the response quantities from zero to t 8 . In 
addition, if it is desired to calculate the response at more than about 120 time 
intervals, it is necessary to make a restart (see Section 6. 1.2). 

For a run from zero, the following indicators are Set: 

RESTRT = 0. 0 

PNCH =1.0 

Here, the condition RESTRT = 0.0 means that it is a start from zero. The 
condition PNCH =1.0 means that at the end of the job certain quantities will 
be punched on cards to be used as data in a future restart. Thus, part of the 
output from this job will be some cards containing the arrays 

TIMX 
ZP(K, L) 

Z2P(K, L) 

Z3P(K, L) 

OMG2(L) 

This punching is done by the executive program 157DR. 

In order to restart the job, the following indicators are set in the 
input data: 


RESTRT =1.0 
PNCH =1.0 

Here, the condition RESTRT =1.0 means that the punched output data of the 
previous job is to be read as input data. The indicator PNCH =1.0 means 
that there will also be punched output at the end of this job. If PNCH = 0. 0, 
no data will be punched and no future restart will be possible. The punched 
cards are put at the end of the data deck. They are read by subroutine 
CRVFIT. 


preceding 


page blank, not bimbo. 


4. 1 INPUT DATA FORMAT 

Data are entered into the program by three subroutines.. The 
executive program 157DR reads the hydrodynamic data, and the DA region 
of the shell data. Subroutine GEOM reads the GDA region of shell data, 
and subroutine CDAFIT reads the CDA region of the shell data. The regions 
DA, GDA, and GDA are read by means of the DECRD subroutine. 
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4. 2 DECRD Subroutine 

The data in regions 
DA, GDA, and CDA is 
read by means of the 
DECRD subroutine, A 
description of the sub- 
routine follows, together 
with a listing in 
FORTRAN IV. 


0 ! 




f — f 

0 


E 


DECRD Decimal Read 


0 6 3 


JLJL3..S. 


2 0 


6 5 


4 . 6. . 4 9 


8 9 6 


E * 0 2 


E 2 


1 1 


E. 2 


2.9 75 


12.3 4 




2 0 0 0 0 0 









2 0 0 0 0 0 


1. Description. When a minus sign is 
encountered in column 1 of a 
DECRD data card, that card will be 
read and then reading will be 
terminated. 

The index of a DECRD card must 
be written to the extreme right of 
the first 12 -column field. 

2. Extent: 78 locations. 

3. Call Statement: 

CALL DECRD (ARRAY) 


The first card will result in informa- 
tion being stored as follows: 


ARR(l) 

ARR(2) 

ARR(3) 

ARR(4) 

ARR(5) 


-0. 7063E 01 
Unchanged 
0. 2435E-00 
0. 2065E 04 
0. 4649E 04 


The - sign in column 1 of the second 
card signals that this is the last card 
to be read under control of this CALL 
DECRD statement. This card has 
been written to illustrate some types 
of errors (or possible errors) in 
writing the data. The information 
will be stored as follows: 


4. 


5 . 


where ARRAY is the name of the 
read array to be read. This argu- 
ment may be subscripted. 

Error indication: If the index field 
is zero or blank, the comment 
"BAD DATA CARD" and the con- 
tents of columns 73-80 of the 
defective card will be printed. The 
job will be terminated. 

Example: Assume a CALL DECRD 
(ARR) statement and the following 
data cards: 


ARR(ll) 

ARR(12) 

ARR(13) 

ARR(14) 

ARR(15) 


0. 7896E 21 (Exponent 
mislocated or 
incomplete. ) 
Unchanged (Treated as 
a blank. ) 

Unchanged (Treated as 
a blank. ) 

0. 2975E 04 
0. 1234E 03 


When no decimal point is written, as 
in the last two items, the data is read 
by the E12. 8 format: the number of 
decimal places is counted from the 
beginning of the exponent field, if any, 
or from the extreme right of the 
field. 
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4. 3 Data Deck Setup 

Data decks should be stacked as follows: 

1. Three title cards (which may be blank, if necessary). 

2. A card with VIN, RHO, WT. 

3. DA, general shell data, read by executive program. 

4. GDA, geometry data, read by GEOM subroutine. 

5. CDA, section properties data, read by CDAFIT subroutine. 

The data in groups 3, 4, and 5 should have a minus sign in column 1 
of the last card. 

The following tables show the nature of the DA, GDA, and CDA decks. 



4.4 Call DECRD (DA) 


DECRD 

Index 

Name 

Description 

1 

EN 

| No, of points along shell meridian 

2 

AO 

Reference length (in. } ) 

3 

HO 

Reference thickness (in. ) | For most cases 

4 

BO 

Reference Young's Modulus (pe\) \ set to 0 

5 

SIGO 

Reference stress (psi) j 

Poisson's ratio 

8 

POI 

11 

SPRL 

Location of spring along meridian 

12 

UK 

Spring value in | direction 

14 

WK 

Spring value in normal direction 

16 

TAU1 

Length of total time interval from zero 

17 

ENTI 

Total no. of time intervals from zero to TAU1 

18 

PI1 

Print interval (will always print last interval) 

25 

MASS 

Mass density lbs. sec^/in^ 

26 

CFE 

Coefficient of viscous damping at each station in 
^direction 

27 

CZ 

Coefficient of viscous damping in normal direction. 

28 

SKFE 

Spring constants of shell under elastic restraint in 
^direction 

29 

SKZ 

Spring constant at each station in normal direction 

30 

SUM 

Fourier summing increment (always -1.) 

33 

TFI 

(Always -1) 

36 

RESTRT 

0. for start from zero, 1. for restart 

37 

PNCH 

0. no future restart; 1. restart cards are punched 

4440 

EMI 

See description of top boundary conditions in Section 4. 6. 

4476 

EMIN 

See description of bottom boundary conditions in 
Section 4. 6. 


Last card must have a - sign in Column 1. 
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4. 5 Boundary Conditions 
4. 5. 1 Top Boundary 

When the boundary conditions on the top boundary are of the following 
kind, a special flag can be used to specify them: 


free: 

(N| 

CD 

£ 

ii 

= ; 

h- M t- 

0) = 

roller: 

(N^ = 

= u e = 

w 

= = 0) 

= 2. 

fixed: 

( v 

: u e = 

w 

= <t»e = o) = 

= 3. 

simply supported: 

(u r 

C 

CD 

H 

w 

= M| = 0) 

= 4. 

complete: 

(u r 

: u 6 " 

A 

F i 

= <(>1=0) 

= 5. 


In these cases, DA(4440) = 1. E10, and DA(4441) is given the value 
1. , 2. , 3. , 4. , or 5. as shown above. Other special boundary conditions 
may also be specified. As an example, the full boundary (which is also given 
above) can be specified as shown in the following data sheets . 
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FORTRAN FIXED 10 DIGIT DECIMAL DATA 

DECK NO. PROGRAMMER DATE ... PAGE 0> 

NUMBER IDENTIFICATION DESCRIPTION DO NOT KEY PUNCH 
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4. 5. 2 Bottom Boundary 

The same selection of boundary conditions is available here as for the 
top boundary. This time, the indicator specifying the free, roller, fixed, 
simply supported, and complete conditions are set as follows: 

DA (4476) = 1. E10 

DA (4477) = 1., 2. , 3., 4. , 5. , 

according to the boundary condition desired. An example of other possible 
boundary conditions is given in the data sheets below. The example here is 
the free boundary (the same as in Section 4. 5. 1). 



FORTRAN FIXED 10 DIGIT DECIMAL DATA 

OCCK NO. PROGRAMMER - OOt <* 
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I •' s 

4/4 CaU DECRD (CDA) 


DECRD 

Index 

Name 

Description 

/ ; 'V/ '** : * _ 

GMI 

Geometry indicator: 



1. « cone - cylinder 



2. = sphere - toroid 



3. a. general discrete point 

4. = arbitrary functions 

2 

EN 

No. of station points 

3 

PFLAG 

Print indicator; 4 0. , prints all data 

4 

RA1 

For GMI = 1. ; radius at station 1 


RC 

For GMI = 2. ; radius of curvature 

5 

AXL 

For GMI = 1.; axial surface length 


ROFF 

For GMI = 2. ; off-set distance to center of curvature 

6 

ANX 

For GMI = 1. ; angle between generator and axis of 



revolution 


prao 

For GMI = 2. ; initial opening angle from vertical axis, 
in degrees. 

7 

PHIN 

For GMI = 2. ; final opening angle from vertical axis, 
in degrees. 

8 

EM 

For GMI = 3. ; number of RI points given 

9-208 

RIPT 

For GMI = 3. ; discrete radii (200 points maximum) 

209-409 

XIPT 

For GMI = 3. ; discrete XI - arc length, (200 points 
maximum) 


The last card must have a - in Column 1. 
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4. 7 Call DECRD (CPA) 

The various tables are set up in this region as follows: 

TAB (1) = No, of stations given along meridian (i, e. , stations 
at which values change), 

TAB (2) s Station No. 1. 

TAB (3) = Parameter value at Station No, 1 
TAB (4) » Next station no. 

TAB ($) * Next parameter value 

I : ; . 

: x-,;-' :■ - ' V '■ " ' , 

I 

Stations and parameter values interlaced. 

The last station must be the Nth station parameter value because 
COOIMA interpolation routine will not extrapolate. 

If +1. 0 x 10^® is placed in TAB (1) the following parameter value is 
constant (uniform over all stations BN) and its value is placed in TAB (2). 


DECRD 

Index 

Name 

Description 

1 

DTB 

Extensional Rigidity 

No* of stations given, if = 1.E10, then a constant 

2 


extensional rigidity is given in 2 
Station No. 1. 

3 

4 

5 

6-41 


if CD A (1) * 1. E10, then this is a constant value of 
extensional rigidity. 

Value Of extensional rigidity between Station 1 and next 
station 

Station No. 2, 

Value of extensional rigidity 

Follow 8 same pattern to DTB (20), value of last rigidity. 



\ 
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DECRD 

Index 

Name 

Description 



Flexural Rigidity 

42 

EKTB 

No. of stations given. if = 1. E10, then a constant flexura 
rigidity is given in 43 

43 


Station No. 1 



if CDA 442) * 1. E10, then this is the constant vajlue of 



flexural rigidity 

44 

-V'w * ; ;. 

Value of . flexural rigidity between station 1 and next 



station 

45 

46 


Station No. 2 

Value Of flexural rigidity 

47-31 

Follows t 

)ame pattern to EKTB (20), value of last ridigity. 
Continue as above for the following quantities: 

83-125 

EITB 

Young's modulus (E) 

124-164 

A.LFTB 

Coefft of thermal expansion (a) 

165-205 

DNATB 

1/2 shell thickness (h/2) 

206-246 

TTB 

Temperature gradient through shell (T) 

247-287 

ENTB 

Membrane thermal load 

288-328 

EMTB 

Bending thermal load 

329-369 

PNTB 

Normal pressure on shell (at reference surface) h| 

370-410 

PFBTB 

Meridional surface pressure (at reference surface) 

452-492 

DZOTB 

Initial displacement in normal direction. 

493-533 

VZOTB 

Initial velocity in normal direction 

534-574 

QZOTB 

Initial acceleration in normal direction 

575-615 

DFOTB 

Initial displacement in £ direction 

616-565 

VFOTB 

Initial velocity in % direction. 

657-691 

QFOTB 

Initial acceleration in £ direction. 


The last card must have a - sign in Column 1. 
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5.1 SAMPLE PROBLEM 

To demonstrate the use of the computer program, and to illustrate 
the format of the input and output data, the sample problem shown in 
Figure 2 has been calculated. 


The problem concerns the vertical impact of a flexible body of 
revolutibn eonsistiog of a shallow spherical shell to which is rigidly attached 
a heavier mast to that their combined weight is 10,000 lbs. The radius of 
curvature of the shell middle surface is 175. 6 ins. , and the opening angle is 
19.53*. The shell extensional and flexural stiffnesses are both set equal to 
3. 33 x 10* lbs/in. , which corresponds to a sandwich shell having 0. 05 in. 
steel facings and 1. 95 in. honeycomb core. Other shell properties are as 
follows: Mass per unit surface area = 9. 7 x 10 lbs. sec. ^/in. 3 ; Poisson's 
ratio = 0. 33, and modulus of elasticity E = 29. 7 x 10^ psi. The initial impact 
velocity is 30 fps. The hydrodynamic loads are computed on the basis of 
the rigid-body theory of Volume 2 of this report, and are then applied as a 
forcing function to the shell of revolution. A full discussion of the numerical 
results obtained is given in Volume 2. Sample data sheets follow. 
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Figure 2. Model of Sample Problem. 
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5,2 Sample Input Data for a start from zero 

In the present problem, we shall use a time interval of . 0001 seconds 
and make a run to 1.0 milliseconds. Because it may, at some future date, 
be desirable to make a restart at 1. 0 ms, (to avoid recalculating the 
response from zero), we shall punch some quantities on cards to make this 
restart possible. They will be part of the input in the future restart. Thus 
we set RESTRT = 0. 0, and PUNCH = 1. 0, 

Sample input sheets are shown below for the start from zero. 
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FORTRAN FIXED 10 DIGIT DECIMAL DATA 

DECK NO. PROGRAMMER DATE ■ PAGE . M 
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PNCH 
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FORTRAN FIXED 10 DIGIT DECIMAL DATA 
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FORTRAN FIXED 10 DIGIT DECIMAL DATA 

DECK NO PROGRAMMER DATE PAGE B» 
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5. 3 Sample Data for the Restart 


In order to restart the problem at 1. 0 ms. , the following input is 
required. Note that the time interval must be the same in all runs. The 
changes in data are indicated by the arrows in the data sheets. Note 
particularly that RESTRT = 1.0 here. The output cards obtained from the 
previous run are placed at the end of the data deck. They are read in 
subroutine CDAFIT. 

The output quantities will be identical to those from the start from 

zero. 
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FORTRAN FIXED 10 DIGIT DECIMAL DATA 

DECK NO PROGRAMMER DATE PAQ P 

NUMBER I IDENTIFICATION I DESCRIPTION DO NOT KEY PUNCH 
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PNCH 




FORTRAN FIXED 10 DIGIT DECIMAL DATA 
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FORTRAN FIXED 10 DIGIT DECIMAL DATA 
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FORTRAN FIXED 10 DIGIT DECIMAL DATA 
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5. 4 Sample Output 


The following pages show the output resulting from the start from 
zero of the sample problem. Letters in circles correspond to descriptions 
below. 

First Title Card 
Second Title Card 
Third Title Card 
No. of Stations EN 
Radius of curvature, ins. 

Offset distance from center of curvature = 0 
Closed apex, therefore PHIO initial opening angle is zero 
PHIN, final opening angle is 19. 53° 

Station number 

R(I), normal distance from shell to axis 
W(THETA) nondimensional curvature in 0 direction 
W(XI) nondimensional curvature in £ direction 
RHOX(I) R(I)/AO = P/AO 
GAMMA (I) p7p 

DTB, the extensional rigidity (constant over shell) 

EKTB, the flexural rigidity (constant over shell) 

EITB, the Young's modulus (constant over shell) 

ALFTB, thermal expansion coefficient (zero) 

DNATB, 1/2 shell thickness (constant over shell) 

TTB, temperature gradient (zero) 
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u 

ENTB, membrane thermal load (zero) 


V 

EMTB, bending thermal load (zero) 


w 

All these other quantities, read by CDR, 

are set to zero 

X 

See Sections 4.4 to 4.7 for descriptions 


EN 

= number of stations 


AO 

= reference length 


HO 

= Reference thickness 


EO 

= Reference Young’s modulus 


SIGO 

* Reference stress level 


ENFO 

= always zero 


ENFL 

= always zero 


POI 

= Poisson’s ratio 


THETA 

= always zero 


PIXI 

= always zero 


SPRL 

= location of spring along meridian 


UK 

= spring value in £ direction 


VK 

= always zero 


WK 

= spring value in normal direction 


EMK 

= always zero 


TAU1 

= total length of time from zero 


ENT1 

= Total no. of time intervals from zero 

to TAU1 

PI1 

= point interval. Here PI1 = 1, and output is pointed at 
end of enc interval 
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TAU2 

= 0 


ENT2 

= 0 


PI2 

= 0 

► (always zero) 

TAU3 

= 0 


ENT 3 

= 0 


PI3 

= 0 


MASS 

= mass density 

CFE 

= coefficient of viscous damping at each station in direction 

CZ 

= coefficient of viscous damping at each station in normal 
direction 

SKFE 

= spring constant of shell under elastic restraining in 
direction 

SUM 

= - 1. 

, always 

EN 1 

= i. , 

open shell; = 2 . , closed shell; set in GEOM 

DEL 

= distance between station points 

BCITP 

= boundary condition indicator - top boundary 

BCIBM 

= boundary condition indicator - bottom boundary 

Y 

Full Tables of 

D 

Extensional rigidity 

EK 

Flexural rigidity 

El 

Young’s modulus 

ALF 

DNA 

Coefficient of thermal expansion (zero here) 
-5“ shell thickness 
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T 


ENT 

EMT 




temperature gradient through shell (zero here) 
membrane thermal load (zero here) 
bending thermal load (zero here) 

Full Tables of PN, PFE, DZO, VZO, AZO, DFO, VFO, 

AFO, all of which were read as zero in CDA. |see Section 4. 7 
for their descriptions J . 

is self-explanatory. It shows the time (t), maximum radius 
of the pressure profile (c), and overall vehicle velocity (V). 
Maximum pressure is given in BC, and depth submerged is 
not computed here. 

is a column showing the station number from the apex (1) to 
the boundary (120). 

is a column of the total pressure acting at each station point. 
Because the maximum radius of the wetted surface is 3. 56 ins. , 
and the distance between each station point (DEL, see fx j) 
is . 503 ins. , the wetted surface only extends to station 87 
Beyond this station, no pressure is applied. Note that the 
maximum pressure always occurs at the edge of the wetted 
surface, in this case at station 8. 

This page contains columns of response output. The columns 
are as follows: 


I 

U(I) 

V(I) 

W(I) 

M(PHI) 

M(THETA) 

M(PHI, 

THETA) 

Q(PHI) 


Station points along shell meridian 

Tangential displacement (ins) of middle surface in £ direction 

Tangential displacement of middle surface in 8 direction. 

Here zero because problem is axially symmetric. 

Transverse displacement (ins. ) of middle surface (positive 
outward). 

Meridional bending moment M^ (in lbs/in. ) 

Circumferential bending moment Mg (in lbs. / in. ) 

Twisting moment M^q. Zero here because problem is axially 
symmetric. 

Shear force Q| (lbs/in. ) 
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Q(THETA) 



I 

N(PHI) 

N(THETA) 

N(PHI, 

THETA) 

SIG(PHI) 

SIG 

(THETA) 

SIG(PHI, 

THETA) 



VEL(U) 

VEL(V) 


VEJL(W) 

ACC(U) 

ACC(V) 

ACC(W) 


Shear force Q 0 . Zero here because problem is axially 
symmetric. 

Additional columns of response. 

Station point along meridian 

Meridional membrane force N| (lbs /in. ) 

Circumferential membrane force N 0 (lbs/ in. ) 

Twisting force N|q. Zero here because problem is axially 
symmetric. 

Stress o-£ (psi) on outer fiber of shell. 

Stress o-Q (psi) on outer fiber of shell. 

Stress on outer fiber of shell. Zero here because problem 
is axisymmetric. 

Additional columns of response. 

Tangential velocity in £ direction (ins. /sec. ) 

Tangential velocity in 0 direction. Zero here because 
problem is axially symmetric 

Transverse velocity (ins/sec. ), positive in outward 
direction. 

Tangential acceleration in £ direction (ins. /sec^) 

Tangential acceleration in 0 direction (ins. /sec^) 

Transverse acceleration (ins. /sec^), positive in outward 
direction. 
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NUMBER OF STATIONS - 120 
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t>. 1 WARNINGS AND RECOMMENDATIONS 


b. 1. 1 Choice of Time Interval 

The proper choice of the appropriate time interval At is important 
lor obtaining good results. If A t is too large, the response will be highly 
damped and inaccurate. On the other hand, if A t is too small, the program 
will take a large amount of time to run. In the sample problem At of 0. 05 ms 
was used, which gave good results. 

o. 1. 2 Number of Iterations 


There is an unidentifiable bug in the problem which makes it necessary 
that a restart be made after about 140 iterations through the shell program. 
Should 160 iterations be exceeded, the program will n blow up. 11 Therefore, 
it is recommended that the job be run in segments of roughly 120 iterations 
between restarts. The method of restarting is explained in Sections 1. 3 and 
5.3. 

0. 1. 3 DECRD 


The subroutine DECRD is in the NAA program library and consequently 
does not appear specifically in the source decks. In installations without 
this program in their library, the subroutine should be inserted in the zero 
link behind the subroutine MMY. 
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